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An experiment to measure neutrons in the upper atmosphere has been 
performed on a balloon flight from Palestine, Texas, at an altitude of 
about 32 km. The experimental arrangement is discussed briefly, and re- 
s-ults of a preliminary analysis of the data for neutrons in the energy 
range 3 to 30 MeV are given. 


The decay of neutrons, produced by the inter- 
action of cosmic rays in the atmosphere and subse- 
quently leaving it (so-called "cosmic-ray albedo 
neutrons"), was suggested initially by Singer 
(ref. 1) as a mechanism to explain the observed 
protons in the Van Allen belt. Hess et al. (ref. 2) 
had measured the spectrum of neutrons in the atmos- 
phere, and these meas\irements were used as the basis 
for subsequent calculations designed to test the 
"Cosmic-Ray Albedo Neutron Decay" (GRAND) theory 
(refs. 3-5). Additional measurements of neutrons 
(refs. 6-9) at high altitudes appeared to confirm 
the spectrum to about 20 MeV since calculated ex- 
trapolations from these measurements agreed with 
the calcxilated leakage flux of Newkirk (ref. 4) 
and Lingenfelter (ref. 5) which use the measure- 
ments of Hess et al. (ref. 2). Note that the ex- 
periments referenced in 6-9 arrive at the "spectrum" 
by assuming a spectral shape and converting it into 
the measured quantity, either pulse-hei^t spectra 
(refs. 6-8) or count rate (ref. 9)« 

Calculations of the proton injection into the 
belt from GRAND neutrons by Dragt et al. (ref. 10) 
using Lingenfelter * s estimate for the neutron 
source, indicated, however, that the GRAND source 
is too small by a factor of 50 for particles with 
energies greater than about 20 MeV. Later measure- 
ments of the proton spectrum in the belt for pro- 
tons above' 10 MeV appear to confirm this discrep- 
ancy (refs. 11-12). While a recent paper by Farley 
et al. (ref. 13) goes far toward resolving the 
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discrepancy, it too uses Lingenfelter ‘ s spectrum 
with the observation that experimented verification 
of this spectrum for energies above 10 MeV is 
lacking. 

It is unfortunately true that the measurements 
by Hess et al. were performed with counters which 
had minimal sensitivity in the energy region from 
rou^ly 10 to 100 MeV. It is in this very region 
that the greatest discrepancies are observed. A 
program -was therefore started in 1966 at the Oak 
Ridge National Laboratory to reduce the uncertainty 
in the neutron spectrum at hi^ altitude in the 
energy region from 2 to 60 MeV, the upper limit 
being chosen because of the availability of detector 
calibration facilities to this energy, and the lower 
limit following from the experimental scheme em- 
ployed. The data-gathering program culminated in 
a balloon fli^t on I7 June 1969 from Palestine, 
Texas. The device floated at an adtitude of 
51,950 ± 200 m (corresponding to about 9 g/cm^ 
residual atmosphere) for about 9-1/2 hours. This 
paper will describe the instrument used and give 
the results of a preliminary analysis of the data 
obtained from this flight. 

The heart of the spectrometer is a liquid 
scintillator (Nuclear Enterprises Inc. NE-215) in 
which neutron interactions eu*e detected by measuring 
the li^t pulses produced by recoil protons. Gamma- 
ray Interactions in the scintillator res\dt in 
pulses which can be separated with pulse- shape 
discrimination (PSD) techniques. Since incident 
charged particles can give rise to p\dses similar 
to those from neutron interactions in the detector, 
the latter is enclosed in a plastic scintillator 
"mantle", of NE-110 (Nuclear Enterprises, Inc., 

San Carlos, California). The output from this 
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’’mantle” is usually put in anticoincidence with the 
output from the liquid scintillator, thus preventing 
any counts due to events which trigger both detec- 
tors. The arrangement used is shown in figure 1. 

The NE-215 is contained in an aluminum can which is 
mounted on an RCA Type -8575 photomultiplier tube. 

The mantle is also viewed by an 8575 tube. At the 
bottom of the assembly can be seen the preampli- 
fier used with the NE-215 detector and the high- 
voltage power supplies for both tubes. The can 
shown on the left is l/l 6 -in. -thick aluminum, used 
as a pressure vessel to maintain the assembly at 
ground pressure to eliminate high-voltage sparking. 
The aluminum can was then surrounded by about 1-1/2 
in. of foam insulation to maintain proper operating 
temperature in the assembly. 



FIGURE 1.— Detector Assembly. 

In order to reduce locally produced background, 
we decided to lower the detector assembly after 
launch about 90 ^ below the gondola on which were 
mounted the electronics and batteries as well as 
misce]J.aneous support equipment. Since it was im- 
practical to retract the detector assembly to the 
gondola prior to landing, the foam insulation was 
in turn contained in a l 6 -gauge sheet- steel can 
designed to take the landing impact. Proof of the 
effectiveness of this arrangement was obtained when 
ye .recovered the equipment and it showed no de- 
terioration in performance. 


A block diagram of the electronics used in 
this experiment is shown in figure 2. Two func- 
tions must be distinguished: pulse-shape dis- 

crimination, and pulse-height analysis. 

As cein be seen from figure 2, the pulse-height 
analysis was conventional. The major item of note 
here is that the discriminators had their thresholds 
adjusted so that the pulse height to be analyzed 
was determined by the threshold of the zero- crossing 
discriminator. 



FIGURE 2.— ^ock Diagram of the Electronics. 


The pulse- shape discrimination was based on 
the difference in the base-line crossing time of 
the amplifier output pulses due to the different 
relative amounts of slow light output from proton 
and electron recoils. The crucial feature of such 
a scheme lies in the proper adjustment and stability 
of the timing involved. The time of an event was 
determined by the input signal to the fast dis- 
criminator connected to the NE-215 photomultiplier, 
and this signal eventually was used as the start 
signal to the time-to-an 5 >litude converter (TAC). 

The TAC stop signal was provided by the zero cross- 
ing discriminator after an arbitrary delay chosen 
so as to put the desired signals into the middle 
of the TAC range. 
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The TAG output was also digitized in a separate 
analog-to- digital converter (ADC); in addition, a 
single channel analyzer used the TAG output to set 
a hit in the telemetry whenever the TAG signal ex- 
ceeded a predetermined level corresponding to p\ilses 
due to incident neutrons. For each event, then, 
the output to the telemetry consisted of the output 
from the two ABG^s and two flags, one to identify 
the PSD ADG and the other to distinguish neutrons 
from gamma rays. A separate telemetry channel 
carried the information from the 8- channel multi- 
scaler shown; yet another telemetry channel, not 
shown, carried housekeeping information such as 
pressures, voltages, and temperatures. 

A preliminary analysis of the data has been 
performed, based on the hardware -set bias for the 
pulse- shape discrimination. A pulse-height spectrum 
of the events so ascribed to incident neutrons is 
shown in figure 3* While a few pulses can be seen 
extending to about channel 250 (corresponding to 
~ 50 MeV incident neutrons), most of the spectrum 
lies below channel 150, or about 50 MeV. 



FIGURE 5.- Pulse-Height Spectrum from Neutron 
Events. 


The analysis performed on these data makes 
use of the calibration of a similar NE-215 detector 
by Verbinski et al. (ref. Ik) and a Spectrum un- 
folding scheme due to Burrus (ref. 15). The cali- 
bration of the detector was extended to about i<-0 
MeV (ref. l6), and a slightly different version of 
the unfolding code (i.e., FERD rather than FERDoR) 
was' employed. The results of this analysis are 
shown in figure k. Also shown on the figure are 
the results of the group at New York University 
(ref. 8) in the 2 to 10 MeV energy range. 

Several points must be made in connection 
with this figure. The first and most important is 
that it is still preliminary; for example, no 
attempt has been made to estimate in any detail the 
required correction for locally produced background 
neutrons. A very crude, zeroth-order calculation 
indicates that such a correction should not amount 
to more than 25^ of the total number of counts ob- 
served, but no detailed calculation to determine 
the spectral shape of this correction has yet been 
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FIGURE 4.— Neutron Energy Spectrum at 9 g/oa® 
Between 95-5° - 100.7® W, and 51.5® - 31.8° N. 
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performed. A second point to note is that the 
analysis procedure results only in the determination 
of a 6 % confidence interval within which the result 
can be expected to lie. It is for this reason that 
no data points are shown but only the limits of 
this confident interval. The width of the interval 
is due not only to the statistical deviations ex- 
pected from the input data but also includes the un- 
certainties associated with the unfolding process 
except those associated with the response matrix. 

For the comparison with the NYU data it must be 
remembered that their data were obtained during a 
period of solar minimum, whereas our data were 
taken during solar maximum. One wo\ild, therefore, 
expect their results to lie above ours, by perhaps 
25 %. 

A comparison with the calculations of Linger- 
felter is best accomplished by comparing our re- 
sults with his in the form of an integral over a 
common energy region. Part of his calculations 
yield O.IO6 neutrons-cm’’^-sec"^-MeV“^ for the 
integrated flux from 5.16 to 10 MeV at 8 g/cm^ 
residual atmosphere for geomagnetic latitude 


(ref. 17). Our preliminary data, integrated from 
3 to 10 MeV for 9 g/cm^ residual atmosphere at 
about 4l° latitude give an upper limit of 0.155 
and a lower limit of 0.142 neutrons/cm“^-sec“^- 
MeV"^. 

Further analysis of the data, as well as an 
atten^t to calculate the neutron spectrum to be 
expected, is in progress. 

We wish to gratefully acknowledge the vital 
contribution to the project made by Mr. H. W. Par- 
ker and his associates from George C. Marshall 
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due Dr. R. B. Mendell of New York University for 
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